QCD Phase Diagram with Imaginary Chemical Potential 

Keitaro Nagata a and Atsushi Nakamura 

Research Institute for Information Science and Education, Hiroshima University, Higashi-Hiroshima 
739-8527 Japan 



Abstract. We report our recent results on the QCD phase diagram obtained from the 
lattice QCD simulation. The location of the phase boundary between hadronic and QGP 
phases in the two-flavor QCD phase diagram is investigated. The imaginary chemical 
potential approach is employed, which is based on Monte Carlo simulations of the QCD 
with imaginary chemical potential and analytic continuation to the real chemical potential 
region. 

1 Introduction 

The QCD phase diagram, which illustrates states of matter formed in terms of the strong interaction 
at a temperature and chemical potential, has been of prime interest in recent physics covering particle 
physics, hadron/nuclear physics and astrophysics. On top of ordinary nuclear matter and hot or dense 
matter such as QGP or compact stars, a very rich structure has been predicted in the QCD phase 
diagram using many phenomenological studies. Thoroughgoing analyses of heavy ion data show that 
we are sweeping finite temperature and density regions. See Ref. |[T|. 

First-principle calculations based on QCD are now highly called. If such calculations would be 
at our hand, their outcomes are also very valuable for many research fields: high energy heavy ion 
collisions, the high density interior of neutron stars and the last stages of the star evolution. Needless 
to say, the inside of nucleus is also a baryon rich environment, and lots of contributions to nuclear 
physics could be expected. 

Because QCD is non-perturbative in most regions of the QCD phase diagram, one is forced to 
use the lattice QCD in order to obtain a quantitative understanding. The lattice QCD is expected to 
provide reliable information on the phase structure based on QCD. Indeed, recently, there have been 
many active quantitative investigations about the finite temperature QCD H2I3H . 

On the other hand, simulations of systems with non-zero quark chemical potential fi have been 
a long challenge for the lattice QCD because of the notorious sign problem. In the lattice QCD, a 
fermionic determinant det/l(yu) is used as a probability in a Monte Carlo method. The introduction 
of non-zero yu makes detA(fj) complex, and therefore leads to the breakdown of the stochastic part of 
the lattice QCD. Despite of the severe sign problem, several approaches have been proposed to study 
the QCD with nonzero yU, where the location of the phase boundary and critical endpoint, and the 
determination of EoS have been extensively investigated. See e.g. H4I5H . 

One idea is to perform simulations in systems with an imaginary chemical potential, where the 
sign problem is absent. The phase diagram in the imaginary chemical potential region is connected 
to that in the real chemical potential region, i.e., ordinary QCD phase diagram, owing to the ana- 
lytic continuation. In addition, HMC algorithms are available without any truncation for the quark 
determinant containing imaginary chemical potential. Numerical costs are relatively small compared 
to approaches requiring the direct evaluation of the quark determinant. Hence, the lattice QCD with 
imaginary chemical potential is one of the standard technique and has been studied by using staggered 
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fermions with two fla vor 11617181 . three flavor 0, four flavor 111011 11121131 in 2-color QCD and finite 
isospin QCD 11141151 . Wilson fermions with two flavor [16| and clover-improved Wilson fermions 
with two flavor IfTTll . In addition, data obtained in such a simulation are also useful the matching of 
phenomenological models such as Polyakov loop extended Nambu-Jona-Lasinio(PNJL) models with 
the lattice QCD II 1912012 11221 . 

Recently, we have studied the two-flavor QCD phase diagram using the lattice QCD simula- 
tions 1181171 . We have used the imaginary chemical potential approach to avoid the sign problem. 
Here, we report our results on the study of the QCD phase diagram. 

2 QCD with Imaginary Chemical Potential 

Let chemical potential complex y. — + //U/(/U e C, Hr,Hi e R). Fermion determinants satisfy a 
relation 

A(rf = y 5 A{-n*)y 5 , (1) 

where A(p.) is a quark matrix. It is straightforward from Eq. (Q]) to obtain (det A(jj))* = detz)(-yU*). This 
implies that detA(p) is complex for a real chemical potential fi = which causes the sign problem. 
On the other hand, one can easily prove detA(fi) is real for a pure imaginary chemical potential fj. - i\xj. 
The sign problem does not occur in this case, and Monte Carlo methods are available. 
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Fig. 1. Schematic figures for the Nf = 2 QCD phase diagram in the (p 2 , T) plane (left) and 
(jj]/T,T) plane (right). A : Pseudo-critical point at ji = 0. B : Critical endpoint. C : Roberge- 
Weiss endpoint. AB : Pseudo-critical line. AC : Extension of the line AB into the imaginary 
chemical potential plane. CD : Roberge- Weiss phase transition line fii/T = n/3. In the right 
panel, larger /ij/T region of the phase diagram is obtained from the RW periodicity. 



Two characteristics of the ft 2 < region are so-called Roberge-Weiss(RW) phase transition and 
Roberge- Weiss(RW) periodicity 11231 . The QCD grand partition function has a periodicity with a period 

2n/N c as 




where k is an integer. Furthermore, Roberge and Weiss showed from a perturbative analysis the exis- 
tence of a first-order phase transition on the line fii/T = n/N c , and from a strong coupling analysis the 
absence of such a transition at low temperatures. Because the RW phase transition occurs at high tem- 
peratures but does not at low temperatures, it may have an endpoint at a temperature Trw on the line 
Hi/T = jt/3. The order of the point has been extensively investigated in Ref. [7|9j and the quark-mass 
dependence is found: first order for small and large quark masses and second order for an intermediate 
quark masses. 
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The left panel of Fig. Q] shows an expected phase diagram, where we employ the (p 2 , T) plane 
containing both the real (/j 2 > 0,p - hr) and imaginary (jx 2 < 0,/i = ipi) regions. Even if /j 2 < 0, 
it is expected quark-gluon-plasma(QGP) and hadronic phases exist at high and low temperatures, re- 
spectively. The two phases are separated by the pseudocritical line for the deconfinement crossover, 
which is an extension from the fi 2 > region. The absolute value of the Polyakov loop is often em- 
ployed to identify confinement/deconfinement phase, although it is not a real order parameter because 
of the crossover nature of the transition. The features in the p 2 < region are well manifested in the 
(jdi/T, T)-phase diagram, see the right panel of Fig.Q] 

3 Result 
3.1 Set Up 

We employ a clover-improved Wilson fermion action of two-flavors and a renormalization-group im- 
proved gauge action. The clover-improved Wilson fermion action is given by 

3 

A{x, y) = 5 xX - [(1 - ydUiWS^ + (1 + 7i )U}(x')6 x ,^ 

i=\ 

- k [e^(l - y 4 )U 4 (x)6 x , tX+i + e^(l + YiWlix^d^] - kC sw 8x,* Yj <7 >° J V 

Here yu is the quark chemical potential in lattice unit, which is introduced to the temporal part of link 
variables. 

In order to scan the phase diagram, simulations were done for more than 150 points on the 
plane in the domain < /// < 0.28800 and 1.79 < J3 < 2.0. All the simulations were performed on a 

x iV, = 8 3 x 4 lattice. The RW phase transition line in the present setup is given by fi; - n/ll ~ 
0.2618. The value of the hopping parameter k were determined for each value of /3 according to a line 
of the constant physics with mp$ /my = 0.8 obtained in Ref. Il25l . 
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Fig. 2. Scatter plots of the Polyakov loop. Left : /3 = 1.80 ( T < T pc ). Right : /? = 1.95 ( 

T > T RW ). 



Scatter plots of the Polyakov loop in the complex plane are shown in Fig. [2] where we choose two 
typical cases /3 = 1 .80 for the hadronic phase and (3 — 1 .95 for the QGP phase. At low temperatures, 
the Polyakov loop is small in magnitude for any fij and continuously changes in a clockwise direction 
as increasing /i/. On the other hand, at high temperatures, the Polyakov loop grows to 0.2 ~ 0.3. It 
stays at the real axis for yn < n/12 and jumps to the left-lower side at fi; - n/12. The difference of the 
Polyakov loop modulus between high and low temperatures shows the deconfinement crossover, which 
is the curve AC in Fig.Q] The observed jump of the Polyakov loop at /i/ = 7r/12 is the Roberge- Weiss 
phase transition, which is the line CD. 
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3.2 Phase transitions and universality 

Now we discuss the properties of the phase transitions in the imaginary chemical potential region of 
the QCD phase diagram. 
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Fig. 3. The /3-dependence of the Polyakov loop modulus L (left) and its susceptibility 
Xl (right) for various /j/. 



The results of the Polyakov loop modulus L and its susceptibility xl are shown in Fig. [3] Those 
behaviors suggest the possibility that the system undergoes the crossover with increasing /3 or tem- 
perature. It is understood from the /ui dependence of the peak position of \l that the pseudo-critical 
temperatures become higher with increasing /// until fij = n/12. This is consistent with the expected 
behavior shown in Fig.Q] 

At fij = 7r/12, the system is on the RW phase transition line at high temperatures. Although the 
behavior of the L andxL are similar, the two state signal is found for fii = tt/12, see Fig. [6] 





MC history MC history 

Fig. 4. Monte Carlo history of the Polyakov loop modulus at ji = (left panels) and ji = n/12 
(right panels). Three values of [} are chosen here. The two state signal is found in the middle of 
the right panels. 
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Fig. 6. Monte Carlo history of the Polyakov loop phase at /ij = n/12 (left panels). The two 
state signal is found in the middle of the left panels. The histogram of (f> at fij = n/12 for various 
fright). 



It is more convenient to consider the Polyakov loop phase <p and its susceptibility x<p m order to 
study the nature of the RW endpoint, which are shown in Fig. [5] For fij = 7r/12, rapidly changes 
near fi = 1.92. It is seen that for fij = n/12 there is one vacuum at low temperatures and are two vacua 
at high temperatures. The histogram and Monte Carlo history of <p at fij - n/12 in Fig. |6]also show 
this behavior. The susceptibility^ shows a divergent-like behavior near /? = 1.92 only for yU/ = 7r/12. 
These behaviors suggest the possibility that the system undergoes the second order phase transition at 
the RW endpoint with increasing temperature. We are now studying a finite size scaling analysis using 
lattice sizes N s = 6, 10. Preliminary result supports that the RW endpoint is second order. 

The /^-dependence of <p and^ are shown in Fig. [7] <p is a smooth function of /i/ at low temperatures 
ifi = 1.80-1.90), while jumps to -27r/3 from at = 7r/12 at a high temperature (J3 = 1.95). 
The system undergoes the first order phase transition at fi[ = 7r/12 at high temperatures. Note L 
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is periodic and <p is anti-periodic, which is caused by the periodicity of the ///-dependence of the 
Polyakov loop GUI . 

Now, the properties of the QCD phase diagram with imaginary chemical potential is summarized 
inFig.|8] 



(1). 




Ml /T 



Fig. 8. Nature of the imaginary chemical potential region of the two-flavor QCD phase diagram 
with an intermediate quark mass obtained from the lattice QCD simulation with the clover- 
improved Wilson fermion. 



3.3 Pseudocritical Line 

Next, we determine the pseudocritical line. Critical values of /? for the deconfinement crossover are 
obtained from the susceptibility of the Polyakov loop modulus for each jjj. We use the result obtained 
by WHOT collaboration E3 to translate p to T. 

We consider quadratic, quartic functions and two types of the Pade approximation 1151141131 : 



1 pc „ X 1 pel 



(3) 
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T vc 1 + £/i(fl//r„ c ) 2 
PI = d w " c; (p ade (I)), (4) 



It \ 

1 pc 

T 
X 1 pc J 



2 , l+di(A 7 /T pc ) 2 



1 + d 2 (fli/T pc ) 2 + diifa/Tpc) 4 



(Pade (II)), 



where fi; = ap./, and /)/ is the imaginary chemical potential in physical unit. Tp C and T pc are pseudo- 
critical temperatures at zero and finite chemical potentials. Note that do(= T pc /T pc (0) at /u - 0) deviates 
from one with 1% because of the disagreement of f3 pc (0) from Ref. [25|. The result is shown in Fig. [9] 
(left panel). The line is almost proportional to the quadratic function at small chemical potentials, and 
shows the rapid increase near /// = 7r/12. We obtain the location of the RW endpoint /3 = 1.927(5), 
which corresponds to T/T pc ~ 1.15. Except for the quadratic function, other three functions are con- 
sistent with obtained data. The quartic function suffers from large errors. The difference between two 
Pade approximations is found near the RW endpoint. 

Obtained pseudo-critical line is analytically continued to /j 2 > region using fi 2 = -fj 2 . Th results 
are shown in the right panel of Fig. [9] The curvature at p./T pc = of a power series of (fi/7rT pc ) 2 is 
?2 = n 2 d 2 - 0.38(12). The present results are slightly smaller than other studies, see e.g. Ref. Il26ll 




4 Summary and outlook 

We have investigated the two-flavor QCD phase diagram using the lattice QCD simulation with the 
imaginary chemical potential approach. The simulation was performed on the imaginary chemical 
potential region. The properties of the imaginary chemical potential region of the QCD phase diagram 
was discussed. We have derived the pseudocritical line, and the results are analytically continued to 
the real chemical potential region. 

The present calculation was performed with the intermediate quark mass and small lattice. The fi- 
nite volume scaling analysis and quark mass-dependence analysis are necessary to confirm the present 
results. In particular, the order of the RW endpoint depends on the mass of the quark. The improvement 
on these points should be done in a future study. The finite volume scaling analysis is under progress. 
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